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A detailed analysis Ms been, carried out of a broadband traveling-wave 
electro-optic light modulator using single crystals of currently available ma- 
terials. The structure analyzed is that first ■proposed by Rigrod and Kaminow 
in which a light beam is reflected back and forth across a microwave trans- 
mission line with the angle between optical and microwave phase velocity 
vectors chosen so that the component of the optical phase velocity vector in the 
direction of microwave propagation is equal to the microwave phase velocity. 
For generality both amplitude and phase modulation are considered in single 
crystals belonging to three different symmetry groups, T d and Did{V d ), 
{linear electro-optic effect) and O h {quadratic effect), for two different orienta- 
tions of the modulating field. 

The inclusion of such practical factors as microwave and optical loss re- 
sults in an optimum design in which the modulator dimensions and operat- 
ing temperature are uniquely determined by the optical and microwave di- 
electric properties of the modulating medium. The main conclusions of the 
analysis are: 

{i) Either cuprous chloride or suitably biased strontium titanate may be 
used to produce 50 per cent linear modulation over a bandwidth of 10 gc 
with less than 5 watts of modulating power, using structures a few centime- 
ters long that have practical manufacturing and operating tolerances. 

{ii) The upper modulation frequency limit is set by a cutoff frequency 
which arises from the finite width of the optical beam. 

{Hi) The presence of appreciable crystalline strain in the modulating 
medium requires that the modulating system consist of a simple phase 
modulator followed by some form of frequency-selective discnminator. Such 
a system requires monochromatic, but not necessarily coherent light. 

I. INTRODUCTION 

Broadband microwave modulation of light requires the use of an ex- 
tended traveling-wave structure. In order to obtain a cumulative inter- 
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action between the microwave and optical traveling waves, some form 
of velocity matching is necessary. In a structure of finite length a prac- 
tical velocity matching criterion is that the component of the optical 
phase velocity vector taken along the microwave wave vector be equal 
to the microwave phase velocity. Several possible arrangements imme- 
diately suggest themselves for realizing this type of interaction. For ex- 
ample, we could use true traveling-wave structures 1 " 4 or iterative struc- 
tures 5 where the optical and modulating signal are adjusted to interact 
intermittently with the proper phase. Broadband modulation further re- 
quires that the modulating structure have little dispersion, which im- 
plies propagation in a TEM or TEM-like mode. Thus, since the dielectric 
constants of usefully transparent modulating media are much less at 
optical frequencies than at microwave frequencies, broadband micro- 
wave modulation of light can be achieved by using some form of optical 
slow-wave structure, of which the simplest is that first proposed by 
Rigrod and Kaminow. 1 In this structure, shown schematically in Fig. 1, 
an optical plane wave zigzags back and forth in an electro-optic medium 
between plane mirrors while the interacting microwave modulating 
signal propagates longitudinally down the structure. The object of this 
paper is to analyze in detail the operating characteristics of this form of 
modulator, to present a systematic and quantitative procedure for de- 
signing modulators that require a minimum of modulating power, and 
finally to present numerical results for a number of presently available 
electro-optic materials. 

The kinds of materials considered here fall into two classes. The first 

of these comprises materials which exhibit a linear electro-optic effect, 

and the second includes materials which possess a quadratic or Kerr 

electro-optic effect. The specific materials considered in the first class are : 

(i) cubic crystals of symmetry T d , such as cuprous chloride (CuCl) 

and zinc sulfide (ZnS), 

(m) KDP and its isomorphs of symmetry D^iVd). 
The second class of materials consists of centrosym metric cubic crystals 
of symmetry O h such as the titanates (SrTiOa and BaTi0 3 ) operated in 
their paraelectric phase. 

All electro-optic modulators using transparent media are basically 
phase modulators in that the modulating fields act directly to change 
the optical phase velocity of the medium. The phase modulation intro- 
duced on a light beam can then be used to transmit intelligence, or if 
desired it may be converted into intensity modulation. This latter process 
may in a sense be viewed as a homodyne process 6 in which we beat two 
phase-modulated normal modes of a birefringent modulating medium 
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Fig. 1 — Basic zigzag modulator configuration (after Rigrod and Kaminow 
Ref. 1). 

(e.g., an ordinary and extraordinary wave), the mixing taking place in a 
polarization analyzer (e.g., a Nicol prism) following the modulator. Most 
amplitude modulators to date have been of this form. 

A systematic comparison of the suitability of the aforementioned ma- 
terials in a zigzag modulator is complicated by the fact that the optimum 
modulator configuration depends upon the form of modulation desired 
(e.g., AM, PM, polarization modulation, etc.) as well as on the demodu- 
lation system used (e.g., simple quantum counter, heterodyne, etc.). The 
comparison is further complicated by the necessity of taking into account 
such practical problems as crystalline strains, losses (both microwave 
and optical) and various operational and manufacturing tolerances. Of 
these effects the existence of strain, particularly random strain, is of 
paramount importance in determining the ultimate form the modulator 
will take. The effect of random strain is to partially destroy the spatial 
coherence of the modulated wave without affecting the temporal co- 
herence. As a result the performance of any modulation-detection system 
which involves the mixing of two optical signals will be degraded, and 
we show in fact that the only system largely unaffected by strain is a 
phase modulator followed by an optical discriminator and quantum 
counter. 

The organization of the remainder of the paper is as follows: In Sec- 
tion II we consider in general the basic phase modulation process in the 
zigzag modulator, obtaining expressions for the time-varying optical 
phase shifts which include the effects of microwave loss and lack of syn- 
chronism between optical and microwave signals. In Section III we use 
these results to obtain expressions for the modulation sensitivity of prac- 
tical phase modulators. In Section IV the same thing is done for ampli- 
tude modulators. The practical implications of such effects as diffraction, 
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strain, optical losses, etc. are discussed in Section V. Finally, in Section 
VI we outline a procedure for selecting the optimum modulator dimen- 
sions and operating temperature so as to obtain the most efficient trans- 
fer of information onto the optical beam. In this same section we con- 
sider the question of modulator bandwidth. 

We conclude with a tabulation of the performance to be expected from 
presently available materials in the zigzag modulator configuration. Our 
results show that it should be possible to build a modulator using either 
the quadratic effect in strontium titanate or the linear effect in cuprous 
chloride that will provide 50 per cent linear modulation over a 10 gc 
bandwidth for modulating powers less than about 5 watts. Two factors 
contribute to the high modulation sensitivity in strontium titanate. 
These are: (i) use of a dc bias results in a large effective microwave elec- 
tro-optic coefficient, and (it) a long optical path is obtained in a small 
volume of material, resulting in a very efficient use of microwave power. 

II. GENERAL THEORY 

The configuration analyzed is shown in Fig. 1. The light is assumed 
to propagate in the yz plane, between mirrors located in the xy plane, 
and to have its electric vector polarized either perpendicular to or parallel 
to the yz plane. The modulating wave propagates in the y direction, along 
the axis of the modulator.* We may distinguish two cases. In the first 
(see Fig. 2a) the electric vector of the modulating field lies in the optical 
plane of incidence (this we call the 0° modulator). In the second (see 
Fig. 2b) the modulating electric field is perpendicular to the optical plane 
of incidence (90° modulator). A TEM modulating wave can be propa- 
gated in the medium with either polarity if the mirrors are made of non- 
conducting dielectric multilayers. 

In the presence of the modulating field and any dc bias, the optical 
index ellipsoid 7 is assumed to take the form shown in Fig. 3(a),f with 
the crystalline axes of the material aligned so that the principal axes of 
the ellipsoid coincide with the coordinate axes.f The significance of the 

* For convenience, we choose to take the xz plane as the transverse plane of the 
modulator and the y axis as the direction of propagation or modulator axis. By 
aligning this Cartesian system with the principal axes of the index ellipsoid of the 
modulating medium, it can be shown that the only effect of changing the direction 
of propagation from, say, the y to the x axis is to cause a sign reversal in the re- 
tardation of the optical wave. _ 

f By this assumption, we exclude cases in which the electro-optic effect takes 
the form of a field-dependent orientation of the index ellipsoid. It can be shown 
that this mode of operation is much less effective than the one in which the ori- 
entation of the index ellipsoid remains fixed, and only the lengths of the axes of 
the ellipsoid change. .... 

t The crystal cut necessary to achieve this alignment for each material is given 
in Appendix A. 
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Fig. 2 — Relative orientation of modulating electric field and optical plane of 
incidence in the two zigzag modulator configurations: (a) 0° modulator; (b) 90° 
modulator. 



index ellipsoid is that it defines the wave vectors and planes of polariza- 
tion of the two normal modes which can propagate through the medium 
in any given direction. A plane perpendicular to the wave vector pass- 
ing through the center of the ellipsoid cuts it in an ellipse whose major 
and minor axes define the directions of the two orthogonal planes of 
polarization. The lengths of these axes are equal to the reciprocals of 
the two indices of refraction. In the configuration shown in Fig. 1, the 
two normal modes have their electric vectors polarized respectively 
parallel and perpendicular to the plane of incidence. The two indices 
of refraction we denote by n\\ and n x ; from Fig. 3(b) it is evident that 

l/n x 2 = 1/n* (la) 
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and 



1 cos 2 Bi , sin 2 6j 



(lb) 



As an optical plane wave traverses the modulating medium, its phase 
velocity varies continuously with time. A convenient way to deduce the 
time dependence of the over-all phase shift, and hence the phase modu- 
lation, is to follow the trajectory of a fixed point on the wave front as it 
moves through the medium. This is similar to the method of character- 
istics used to treat nonlinear and parametric wave propagation problems 
and is valid whenever the phase of the wave or the eikonal of geometrical 
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Fig. 3 — (a) Orientation of index ellipsoid with respect to coordinates of modu- 
lator; (b) projection of index ellipsoid on optical plane of incidence. 
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optics is largo. When this condition is satisfied the wave propagates 
approximately as a plane wave. This can be demonstrated by expand- 
ing the eikonal in a power series over small space regions and time inter- 
vals. 8 It is then correct to say that the time taken for a point on the 
wavefront to advance a distance y along the modulator axis is given by 

(I - , ) = J™ (2) 

c sin di 

where c is the velocity of light in vacuo, n a is the refractive index in the 
absence of any electric field, / is the starting time at y = 0, and t is 
the angle between the light ray and the z axis. This simple expression 
replaces a complicated integral equation and is correct to first order in 
the fractional change An/n in the index of refraction, and hence is valid 
only when | An/n | « 1, as is usually the case. 

The value of the wave phase at any point in space is given by the 
integral 

where io/2ir is the optical frequency, and f is measured in the ray direc- 
tion. This equation can be written in terms of the component of the 
phase velocity vector in the ^-direction* as 

sin di Jq Vp(y,t) Xo sin 0, <>o 

where X is the free-space optical wavelength and S the modulator length. 
Equation (3) applies separately to both normal modes of the electro- 
optic medium when the appropriate index of refraction n is used. How 
n depends on the applied electric field for the various materials and 
configurations of interest is derived in Appendix A. 

By substituting (2) into (3) and utilizing the expressions derived in 
Appendix A for n(y,t), one can calculate the time dependence of the 
wave phase correct to first order in the fractional change An/n. As long 
as the modulator length is bounded, as it must be if diffraction effects, 
etc., are to be minimized (see Section V), this calculation should give 
accurate results since | An/n | « 1. 

In the remainder of this section we make use of (2) and (3) to obtain 
expressions for the wave phase as a function of time for the two normal 
modes of the modulating medium. These expressions describe the basic 

* Note that when dispersion can be neglected the y component of the optical 
phase velocity vector is just the group velocity of the optical energy flow along 
the modulator axis. 
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modulation process within the medium. In subsequent sections we shall 
go on to describe the practical considerations involved in the actual 
realization of phase and amplitude modulators. 

2.1 Linear Electro-Optic Effect in Cubic Materials 

2.1.1 0° Modulator 

For this case the two modes have refractive indices given by ( 1 ) and 
(104) 

- = - (l + n.V«*.)* (4) 

Wj. no 

- = -(1 - n, 2 r*E z cos 2 •«)* (5) 

w n no 

where r« is the linear electro-optic coefficient (see Appendix A). For an 
attenuated modulating wave of the form 

E,(y,t) = E m e~ amV cos (uj - /3„,/y) (6) 

the phase of the wave polarized perpendicular to the plane of incidence 
is approximately given by 

<P x (to) 

w 2* _n«_ [* n_ ^ ni E m e~ amV cos [*Uo + (jB. - /O//]} % (?) 
Xo sin 0i Jo 

where for convenience we have introduced the notation 



c sin 0, 



- ft . (8) 



The approximation holds when \n 2 r»E m \ « 1. Equation (7) can be 
integrated directly, giving 

rxw Xo sin 0, Xo sin Oi 

X ^ aJ + fjt - ) cos («J + *) 

where /3 is defined by 

j8 = 0. - /3™ (10) 
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and \p is a time-independent phase factor. In the absence of microwave 
losses this reduces to the result obtained earlier by Kaminow and Liu. 2 
It is evident that the phase modulation described by (9) is maximized 
by the choice /3 = /3 e — /3 m = 0, corresponding to the matching of the 
^-components of the microwave and optical phase velocity vectors. 
Thus, it is clear that microwave losses do not affect the basic velocity 
matching condition. Under this condition (9) reduces to 

{t) = 2tt noS {l _ L n 2 riiEmf(amS/2) cos Uat] (u) 

Xo sin 6i 

where for the sake of brevity we have defined the function f(a m S/2) by 

-(-¥)P^>'(¥)- <» 

For the other wave, whose plane of polarization lies in the plane of in- 
cidence, the corresponding expression is 



„,(0 . ^r-JkjJ [i + i n * riiEm cos 2 dif /W>\ cog u l (13 
Ao sin at [_ \ z / j 



2.1.2 90° Modulator 



The analysis for the 90° modulator is formally identical, except that 
(105) is used for the principal indices of refraction, in place of (104), so 
that 

l/n x = l/n (14a) 

and 

— = - (1 - noruE t cos 26 i )\ (14b) 

n n n 

The results for a modulating wave of the form given by (6) are 

ipji) = z : — - = constant (15) 

Xo sin di 

and 

Vl(t) - ^ -?4- T 1 + WnxE m cos 2Bif fa£) cos w m t] . (16) 
Ao sin ui\_ \ z / j 

2.2 Quadratic Electro-Optic Effect in Cubic and /or Amorphous Materials 

The analysis for the quadratic effect in the presence of a large dc bias 
is very similar to that for the linear effect. The situation is complicated 
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somewhat by the presence of two independent electro-optic coefficients 
in the case of the quadratic effect in cubic crystals (see Appendix A). 
We consider the two configurations discussed before, in which the modu- 
lating field lies parallel (0° modulator) or perpendicular (90° modulator) 
to the optical plane of incidence. 

2.2.1 0° Modulator 

For this case the relevant refractive indices are [refer to (110)] 

- - - (1+ a. 1 **. 1 ) 1 (17) 



no 



and 



-1 = I [i + n > pi E> (sin 2 6i + ^ cos 2 d)l (18) 

w n n L \ Pi /J 

where p x and p 2 are quadratic electro-optic coefficients (see Appendix A). 
In order to achieve approximately linear modulation, as well as to en- 
hance the modulation sensitivity, the applied field should consist of a 
large static bias and a smaller microwave field, so that 

E t (y,t) = E + E m e~" mV cos (coj - p m y). (19) 

Then to first order in E m /E Q 

E z 2 tt E 2 + 2E E m e~ amV cos (o>J - p m y) 

and the expressions for l/n ± and 1/ny become linear in E m . The inte- 
grations needed to obtain <p ± and <p\\ are similar to those in the preceding 
section. Again the maximum ac phase shift occurs with velocity match- 
ing, for which case we obtain 

„ A (0 = 2 JL^LU- ± no 2 P ,E [#o + 2E m f (°^f) cos «„*]} (20) 
and 



^ ( t ) = 2 JL ** U - i 2 n 2 Pl E (sin 2 Bi + * cos 2 d) 
Xo sin 6i { \ Pi / 



X [#o + 2/iv(^)cosauJj . 



(21) 



In deriving (20) and (21) we note that because the dc bias makes the 
material birefringent even in the absence of a microwave modulating 
field, the ordinary (J.) and extraordinary (||) modes cannot, strictly 
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speaking, be simultaneously matched to the modulating wave. We as- 
sume, however, that the bias is small ( | n 2 pE 2 | « 1 ) so that the dis- 
tinction between ordinary and extraordinary phase velocities becomes 
unimportant. For known materials this condition is met for fields as 
large as 10 5 v/cm — i.e., for any field that the material can sustain 
without dielectric breakdown. 

2.2.2 90° Modulator 

The analysis here is similar to the 0° case, except that the indices of 
refraction are simply [refer to (111)] 

- = - (1 + noW?, 2 )* (22a) 
n ± n 

and 

- = -(1 + noW?* 2 ) 1 (22b) 
n o wo 

so that the resulting wave phases are 
, x (0 = *E »£- |l - Wp^o [tf. + 2E m f (^fj cos w j| (23) 

and 
<P, (0 = j£ ^ (l - WW?o [^o + 22T m / (^) cos ou]} . (24) 

2.3 Linear Electro-Optic Effect in Uniaxial Materials 

The results here are similar to those obtained earlier for the effect 
in cubic materials. 

2.3.1 0° Modulator 

Equations (1) and (113) apply here, so that 

— = - (1 +no*r»#.)* 

1 ~ l A 2 „ 2o 2An . 2fl Y 

— « — I 1 — n r 6 3#* cos t - — sin 0,- ) 

nn n \ n / 

where we have set An = n e — n , which we assume to be much less 
than unity (n e refers to the extraordinary refractive index). The wave 
phases for the two modes are then given by 
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t) = 2jr ^ojS r i no V 6A/ /o^S\ cog H (26) 

Xo sin 0,- L \ 2 / J 



and 

,,v 27T WoS 

" (0 = x75nT < 



X 1 — - — sin 2 6i — %n<?r«&E m cos 2 t f (^— ) cos u m t . 



2.3.2 90° Modulator 

From (1) and (114) we obtain 

\/n x = 1/n 
and 

1 1 



so that the wave phases are 



and 



(27) 



= — (1 - noVeatf* cos 2 0<)* (28) 

nil n 



<p,(t) = :— — = constant (29) 

Xo sin 6i 



vi (0 = v ^4" T 1 - WrtoB. cos 20,-/ (^) cos uj] . (30) 
Ao sin ui]_ \ £i / J 



III. PHASE MODULATORS 

The preceding results are immediately applicable to the phase modu- 
lator. The maximum modulation sensitivity is obtained when the beam 
is linearly polarized in the direction ( _]_ or || ) showing the greatest effect, 
although unpolarized light can also be used. 

For all the cases discussed in the previous section, the ac component 
of the wave phase can be written in the form 



ASgMf (-: 



*..(*) = 2j P~ ASgMf (?g) E m cos o> m t, (31) 



where A is proportional to electro-optic coefficient and g(6i) depends 
on the angle of incidence. The quantity of most immediate interest is 
the modulation sensitivity, which we may define as the ratio of the 
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peak phase shift to the square root of the power input to the modulator.* 
To make the problem definite, we consider a parallel-plate line of spac- 
ing h and width w, filled with a medium of microwave dielectric constant 
k. If we ignore fringing fields,! the power flow is given in terms of the 
electric field by 

P = \y^. E Jhw, (32) 

where 77 = 377 ohms is the impedance of free space. The modulation 
sensitivity is then 



M 

Vp 



2irti 3 



^ ( *>'(¥)l/=S; ™ 



where M is the peak ac phase shift. Since the angle of incidence is deter- 
mined in terms of material parameters by the velocity matching condi- 
tion /3«. = ft m , which reduces to sin 0,- = Wo/vV for the present geometry, 
(33) can also be written as 



M 



VP 



2 ¥^&)^)Vvk» <*> 



where pino/yAc') is a function of (rio/vV') which replaces g(d { ). Values 
of A, g{6i) and p(rio/vV) for the cases of interest here are given in 
Table I. In this table, the linear electro-optic coefficient which we write 
simply as r,y is to be taken to be ni in the case of cubic materials, and 
r fl3 for uniaxial crystals. 

Numerical results for the power required in practical cases are given 
in Section VII. 

IV. AMPLITUDE MODULATORS 

In this section we treat the problem of obtaining broadband ampli- 
tude-modulated light by suitably combining the two simultaneously 
phase-modulated normal modes of a transparent electro-optic medium. 
In the ensuing treatment we will consider only the intensity variations 
in the modulated beam. In so doing we are ignoring the phase informa- 
tion that would be recovered, for example, by an optical heterodyne re- 

* The advantage of defining a modulation sensitivity in this way is that it is a 
function of the microwave-optical properties of the electro-optic modulating me- 
dium and is expressed in terms of power which is easily measured. 

t The effects of fringing fields and the problems associated with non-TEM mode 
propagation in this structure are considered in Section 5.1. 
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Table I — Functions Appearing in (33) and (34) for 
Various Phase Modulator Configurations 



Modulator 


Mode 


A 


S(0i) 


/>(»o/\/«') 


Linear, 0° 


_L 
II 


|r., 


1/sin 0, 

COS 2 di 

sin di 


yV / _ ftf\ 

no \ «' / 


Linear, 90° 


-L 



If... 


cos 20, 
sin di 


no \ »c' / 


Quadratic, 0° 


_L 


pzEo 
PiEo 


1/sin 9i 

P2 

sin 2 Oi + - cos 2 di 
Pi 


VW"o 
k' pi \ K ) 




sin 0, 


n 

V7' 


Quadratic, 90° 


1 


pil'J a 
P2E0 


1/sin 6i 


Vx'/rio 



ceiver, and in effect are assuming that the receiver is a simple quantum 
counter. Thus we look for modulation of the form 



I(t) = J (l + M COSWmO 



(35) 



where Iq is the average light intensity and M is the modulation index. 

Cubic materials are of particular interest as amplitude modulators 
since they are isotropic in the absence of an applied field and become 
uniaxial only upon application of a field. This is in contrast to what 
takes place in materials with lower degrees of symmetry where naturally 
occurring birefringence exists. The significance of starting with an iso- 
tropic crystal, as opposed to a uniaxial or even biaxial crystal, is that in 
the former no static retardation is accumulated between the two normal 
modes. Hence any ac difference in phase produced between these modes 
can easily be converted into useful intensity modulation. The presence 
of a static contribution to the differential phase shift can cause the in- 
duced ac birefringence to be masked out due to the angular spread of 
the beam, constructional uncertainties or finite source linewidth. (See 
Section 4.3.) 

There are two simple ways in which a time-varying phase shift be- 



LIGHT MODULATORS 



2G35 



tween two orthogonally polarized light beams may be converted into 
approximately linear intensity modulation. One of these is most useful 
when there is no dc retardation between the two beams, and the other 
when there is a controllable dc retardation, as in the biased quadratic 
effect modulator. 



4.1 Amplitude Modulation Using the Linear Electro-Optic Effect in Cubic 
Materials 

The first method of producing amplitude modulation is shown in 
Fig. 4. Circularly polarized light is produced by passing the linearly 
polarized incident beam through a quarter-wave plate in the manner 
shown. 

The light emerging from the quarter-wave plate may be decomposed 
into two orthogonal normal modes, one aligned with the fast or f axis 
of the quarter-wave plate and the other aligned with the slow or £ axis 
(see Fig. 4). Both normal modes are of equal amplitude but separated 
in phase by 90°. In order to obtain the desired modulation, the plane 
of transmission of the analyzer must make an angle of 45° with the plane 
of incidence. 

For the arrangement given in Fig. 4 one can easily show that the ratio 
of transmitted to incident light intensity is given by 



It/h = |[1 + sin f> A - <e|)]. 



(36) 




POLARIZER 




A FAST(f) 
^|45 oX AXIS PR| SM 



QUARTER-WAVE 
PLATE 



MODU 



PRISM 



ATOR 



Fig. 4 — Arrangement for obtaining linear amplitude modulation for cubic 
crystals exhibiting a linear electro-optic effect. This arrangement can also be used 
With uniaxial crystals if the quarter-wave plate is replaced by a suitable optical 
compensator. 
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For the modulators using the linear electro-optic effect in cubic crystals, 
(11), (13), (15), and (16) yield 



<p± 



<Pu = 



Xo 



EnJSgMf 



(¥) 



COS U) m t 



(37) 



where 



g(6i) = 



1 + cos 2 Oj 



sui 6 
cos 20,- 



sin Oi 



(0° modulator) 
(90° modulator). 



(38) 



To obtain a specific expression for the modulation sensitivity, which we 
define as | M/\/P | as before, we may again consider the simple parallel 
plate geometry discussed in Section III. It is evident from (36) and (37) 
that the modulated signal will in general consist of an infinite number of 
sidebands whose amplitudes are determined by Bessel functions. For 
simplicity we limit ourselves to the case of small modulation index (viz., 
I <P± — <P\\ I < t/6) so that only the first sideband is important. We 
then obtain for the modulation sensitivity 



M 



Vp 



_ wVv/2 S 



V\ Vh 



w 



*'&)'(*£)• (39) 



with 



te) 



2 - n„7/c' 



Wo/Vk' 
1 - 2ti 2 /k' 



n 



/Vk' 



(0° modulator) 
(90° modulator). 



(40) 



In this expression we have again used the velocity matching condition 
sin di = uo/x/k' to replace g($ t ) by p(n /\/ii'), and have replaced the 
electro-optic coefficient with the more generally useful half-wave retar- 
dation voltage, V\ , defined by* 



V, = 



Xo 



2no 3 r 4 i 



(41) 



Once again we defer computing actual numbers until Section VII. 



* Note that when applied to the transverse electro-optic effect V\ need not be 
the actual voltage applied to the crystal to obtain ir phase retardation, which de- 
pends on the geometry of the crystal, but iB to be viewed as a convenient collection 
of constants. 
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4.2 Amplitude Modulation Using the Quadratic Electro-Optic Effect 

Use of the quadratic effect in conjunction with a dc bias to provide 
amplitude modulation differs from the linear case just considered in 
that there is now a bias-field dependent static phase difference produced 
between the ordinary and extraordinary rays. Approximately linear 
intensity modulation can be obtained in this case by placing the active 
medium between crossed polarizer and analyzer, as shown in Fig. 5, 
provided that the dc bias is adjusted so that operation occurs about the 
point of half-transmission. For the arrangement of Fig. 5, the ratio of 
transmitted to incident light intensity is given by 

UIU = \\\ + cos { Vx - *)]. (42) 

According to (20), (21), (23), and (24), the phase difference for the 
quadratic effect can be written as 



<p± — <P\ 



irtlo 



(pi - ptWgidi) \e + 2E m f fe^j cos w m t\ (43) 



where 



g(o<) = 



sin 0, (0° modulator) 

1/sin Bi (90° modulator). 



(44) 



Equation (43) is in the form <p ± — py = 7d + Tao ; hence (42) can be 
written : 



It/h = Ml + cos Tdc cos 7 ac — sin 7 do sin ? ao ]. 



(45) 




POLARIZER 



MODULATOR 
CRYSTAL 



Fig. 5 — Arrangement for obtaining linear amplitude modulation in materials 
exhibiting a quadratic electro-optic effect. 
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For 7ac < x/6 (small modulation) this becomes, approximately 

Il/h ~ §(1 - Tac shiTdc). 



(46) 



Maximum modulation sensitivity results when the operating point is 
chosen so that* 

TTWo (pi — P2)E Sg(di) 



Tac = 



- (2w+l)£, 



m = 0,1,2,.--. (47) 



The point that should be noted here is that for the linear Kerr effect 
modulator there are modes of operation and that these modes, which 
are designated by the index m, constitute an independent operating 
variable. 

For the operating points defined by (47), the modulation sensitivity 
for the parallel-plate geometry can be written in the form 



M 



Vp 



-^AM^Wt)- (48) 



where 



Mo 



no 



no 



(0° modulator) 



(90° modulator). 



(49) 



In deriving this equation we have, as before, made use of the synchronism 
condition sin 0* = Wo/Vk 7 , and have replaced the electro-optic coefficient 
(pi — Pi) by the more familiar Kerr constant 



K = 



n (pi — pi) 
2X 



(50) 



It should be noted that the dc field is not an independent variable, but, 
because it determines the operating point [see (47)], is limited to the 
discrete values 



*-i 



(2m + 1) 



4/ KSpC 



no 



2 n^y 



(51) 



In practice E would be made as large as the dielectric strength of the 
material would permit, subject to (50). We should point out that the 

* It is obvious that the arrangement shown in Fig. 4 could also be utilized here 
if desired; the only difference would be that the quiescent operating point would 
be chosen to make the dc part of <p , — <p u in (43) an odd multiple of w. 
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effects of electrostriction at high dc fields will undoubtedly result in 
electrostrictively induced birefringence. To the extent that the electro- 
striction is uniform, it is equivalent to a change in the Kerr constant, 
and can be compensated for by timing E . On the other hand, nonuni- 
form electrostriction can degrade the modulator performance in the 
same way as nonuniform strain. (See Section 5.6.) 

4.3 Amplitude Modulation Using the Transverse Linear Electro-Optic 
Effect in Uniaxial Crystals 

Use of the zigzag configuration to produce amplitude modulation 
with noncubic materials is complicated somewhat by the presence of a 
large static phase retardation due to the natural birefringence of the 
material. When the modulator is used with well-collimated monochro- 
matic light, however, the effects of the natural birefringence can in 
principle be eliminated by means of an optical compensator, thereby 
making it worthwhile to consider this case as well. 

The physical arrangement used in conjunction with the linear effect 
in cubic crystals, and shown in Fig. 4, is also applicable here if we re- 
place the quarter-wave plate by a suitable compensator. The compen- 
sator allows the phase of one optical mode of the modulator to be varied 
relative to that of the other. 

For this arrangement the ratio of transmitted to incident intensity 
is given by 

h/h = |[1 + cos (<p ± - n + *>„)] (52) 

where <p c is the relative phase shift introduced by the compensator. The 
phase retardation in the modulator is given by (26) to (30), and has 
the form 

(f>x - n = - ^ 3 E m Sg(O t )f fef) cos («J) + » (53) 

where g(6i) is given by (38) of Section 4.2 and <p„ is the static retarda- 
tion, given by 



<Ps = 



sin 6i (0° modulator) 

Xo 

2tt — -J— ( 90° modulator) . 



(54) 



Xo sin 8i 

It is evident from (52) and (53) that maximum linear amplitude modu- 
lation is obtained when the compensator phase shift is adjusted so that 

<p c + <p s = (2k + 1)tt/2, k = 0, 1, 2, 3, • • •, (55) 
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for which case (52) becomes 

It/h = Ml - ("l) (i " +1) sm (<p ± - n ) ac ]. (56) 

With this adjustment of <p c , uniaxial and cubic materials are optically 
equivalent as far as modulator operation is concerned [compare (53) 
and (56) with (36) and (37)], so that (39) to (41) for the modulation 
sensitivity can be applied here as well by changing r 4 i to n 3 . 

The question arises, of course, as to the practicality of compensating 
for the natural birefringence over a finite aperture; if the modulator is 
to function properly, the static retardation <p 8 must be the same, within 
a fraction of a radian, for all rays comprising the beam. Three factors 
are involved here: 

(i) monochromaticity of the optical source, 
(it) dimensional tolerances, and 

(Hi) angular spread in the incident beam. 
To illustrate the problem we consider the 0° modulator in KDP. In this 
configuration the static retardation is given by 

<p a = 4tt — — sin 0,- . (57) 

Ao 

If we neglect dispersion in An, a range of wavelengths AX in the source 
leads to a spread in <p 8 of 



A<p 8 I tt <p s 



AXn 



Xo 
while a variation AS in the modulator length results in a spread 



(58) 



A<p a tt cp 8 — . (59) 

Under typical circumstances (An = 0.04, 5 = 4 cm, 0» = 19°, X = 0.6 n; 
see Section VII) <p„ is of order 10 4 radians, so that the requirement 
A<p s < 1 implies that the source monochromaticity and modulator di- 
mensions be maintained within a few parts in 10 5 . Although the fore- 
going tolerances are tight, they are attainable. 

A more fundamental problem concerns the angular spread in the in- 
cident beam as it traverses the modulating medium. For example, let us 
suppose that the beam divergence can be characterized by a spread A0,- 
in the angle of incidence. This will produce a variation in <p 8 given by 

A<p 8 « 4tt ^ cos BiAdi (60) 

Xo 
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so that the inequality Ay„ < 1 requires that 

M t < —- £ . (61) 

4irA)iS cos 8i 

For the example presently under consideration, (61) gives A0; < 3.2 X 
10 -B radian = 6.5 sec of arc. Whether or not this requirement can be met 
depends on the details of the optical system. If, for example, the modula- 
tor were operated in the far field of an optical maser of output aperture 
/) which is smaller than the modulator aperture, diffraction alone would 
produce a spread A0, ~ 1.22a /-D, so that the inequality could be satisfied 
only if the beam diameter were greater than 2.3 cm, which is clearly 
unreasonable. On the other hand, if the modulator were placed in the 
near field of the optical maser, (which extends for a distance D /X , 
which is of order 1 m for a 1-rara diameter beam), diffraction could be 
negligible, and strain refraction would then set a lower limit to the mini- 
mum realizable effective A0,- , which would in all likelihood still be con- 
siderably in excess of the limit set by (61). 

In view of the foregoing practical considerations it would seem un- 
reasonable to attempt the construction of an amplitude modulator in 
the zigzag configuration using a material such as KDP. Moreover, as 
will be shown in connection with the. KDP phase modulator, this ma- 
terial typically requires at least an order of magnitude more power than 
other materials under consideration. 

4.4 Summary of Results 

The results for the various forms of amplitude modulator may be sum- 
marized in a single equation, of the same form as that obtained for the 
phase modulators 

Values of A and p(/io/vV) are given in Table II. 

V. FACTORS AFFECTI.N'G PRACTICAL MODULATOR PERFORMANCE 

It is appropriate before discussing the actual design and performance 
of zigzag modulators to consider some of the practical factors which 
affect the operation of the device. In particular we consider the effects of 
(i) fringing fields 

(ii) deviations from the correct angle of incidence 
(in) diffraction 
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Table II — Functions Appearing in (62) for Various 
Amplitude Modulator Configurations 



Modulator 


A 


g(8i) 


p(no/Vx') 


Linear, 0° 


fr« 


1 + cos 2 0i 
sin 0i 


2 - «„«/«' 
n /vV 


Linear, 90° 


\Tii 


cos 20, 
sin 0i 


1 - 2no I A / 
«o/VV 


Quadratic, 0° 


VPl — P2)Eo 


sin Bi 


no/VS 


Quadratic, 90° 


(pi — P2)Eo 


1/sin 0i 


Vic' /no 



(iv) optical losses 

(u) lack of parallelism in the reflecting surfaces, and 
(vi) strain. 

5.1 Effect of Fringing Fields 

The calculations made thus far are strictly correct only in the limit of 
a uniform microwave field across the active cross section of the modula- 
tor. Whether this condition is satisfied or not depends upon the ratio of 
the width of the electro-optic modulating medium, w, to the width, H, 
of the parallel plate conductor (see Fig. 6). This ratio, given by 

£ = w/H (63) 

determines the filling factor of the dielectric loaded parallel plate guide. 
It can be shown from the analysis of Kaminow and Liu 2 that when £ = 1 
the field distribution in the modulating medium is uniform and propa- 
gates as a TEM wave. This result is, however, valid only when open- 
circuit boundary conditions prevail at the ends of the waveguide cross 
section. 2 These boundary conditions will apply physically only to the 
extent that fringing fields can be neglected. This condition can be ex- 
pected to be satisfied when the change in transverse impedance in going 
from the dielectric region to air is large, or alternatively, when the energy 
stored in the dielectric is large compared to the energy stored in the 
fringing fields. As long as the dielectric constant of the electro-optic 
medium is in excess of 10, fringing fields should be negligible and there- 
fore have relatively little effect on the foregoing calculations. 

It is important to note that it is generally difficult to propagate TEM 
or TEM-like waves over a broad band of frequencies in partially loaded 
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two-conductor transmission lines of the form shown in Fig. 6. This is 
because dielectric waveguide modes can be supported in such structures 
at high frequencies. The relative importance of the dielectric modes de- 
pends on the filling factor £. For filling factors less than unity, three 
propagating regions can be distinguished. 2 The first occurs at low fre- 
quencies where 

X,„ » Vk'w. (64) 

In this region the field distribution in the dielectric region is uniform and 
TEM-like mode propagation results. If, on the other hand, we have 



X™ « vV 



w 



(05) 



then, as shown in Ref. 2, the microwave energy resides primarily within 
the dielectric with a cosinusoidal field distribution, and once again 
propagates as a TEM-like wave. For those cases where X„, ~ v k'w we 
find that a symmetrical TEio-like waveguide mode propagates with the 
dispersion characteristic of such a mode. 

What is of primary significance here is that for values of the filling 
factor between 

0.5 < £ ^ 1 

the dispersion in the microwave phase velocity is small (from the calcula- 
tions given in Ref. 2 a 30-40 per cent variation in phase velocity results 
for £ = 0.5). Moreover, as the filling factor approaches unity, not only 
does the dispersion vanish but the field distribution also becomes uni- 
form. Based upon these considerations it is evident that a parallel plate 
guide completely filled with a moderately high dielectric constant elec- 



H 



r 
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Fig. 6 — Representative microwave electric field distribution in the transverse 
plane of the modulating medium. 
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tro-optic material is a suitable practical broadband microwave guiding 
structure. 

5.2 Sensitivity of the Modulator Performance to Small Changes in the Angle 
of Incidence 

The velocity matching condition requires that the angle of incidence 
be adjusted so that sin 0, = n /vV. To determine the effect of small 
deviations from this value, we return to the general expressions for the 
basic phase modulation processes [(9) and following] and observe that 
in the lossless case (a m = 0) the modulation sensitivity is always of the 
form 



M 



:-0S 



Vp 



sin 



-= Cgiit) -^- (66) 

T 



where g($ t ) is a slowly varying function of 0, , and the parameter /3 is 
given by [see (8) and (10)] 

p Pm \ sin 6i ] 

Neglecting any variation of g(di), we find that a variation in 0, of 

\Ut\*0A%$(l-?£f (67) 

from the value for perfect velocity matching will reduce the modulation 
sensitivity by a factor of 2. 

5.3 Diffraction Effects 

Diffraction can limit the performance of the zigzag modulator by pro- 
ducing a divergence in the beam as it passes through the modulator. This 
is equivalent to a spread in the angle of incidence which, for satisfactory 
modulator operation, must be less than the value of A0, given by (67). 
Whether or not this requirement can be met depends primarily on 
whether the modulator is operating in the far or near field of the optical 
source. The near field of an aperture of diameter D b illuminated by a 
coherent plane wave extends a distance of order Db Ao • In this region, 
which is about 1 m for a 1-rara diameter beam, there is little spread. At 
distances appreciably greater than this, the beam propagates according 
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to the simple theory of Fraunhofer diffraction, with the main lobe having 
an angular width 

"~ L22Xo (68) 



Thus, in the worst possible case diffraction could result in an effective 
spread in the angle of incidence of order 1.22\ /D b , and hence limit the 
allowable value of the ratio of modulator length to beamwidth, accord- 
ing to (67), to 

fljij^A *V (69) 

Dd 2 Ao K \ K / 

The aperture D b is further restricted, in terms of the geometry of the 
modulator, by the requirement that the beam remain in the modulating 
medium after the first bounce. From Fig. 1 it is evident that this requires 
that 

D b g 2d sin 0, , (70) 

where d = h in the 0° modulator, and d = w in the 90° modulator. 
Equations (69) and (70) together require that 

o ^ Xm no 1 (7\) 

d = \o k'I (1 -no 2 A') r 

5.4 Effect of Optical Losses 

We shall distinguish between two kinds of optical loss. The first arises 
from scattering due to crystalline inhomogeneities and absorption within 
the bulk material, and is described by an absorption coefficient a s . The 
second is the reflection loss at the mirrors, which are assumed to have a 
reflectivity R. The optical transmissivity can then be written as 

T/h = exp [- asS esc 0, + (S/d) cot 0, In R] (72) 

when (S/d) cot 0,, the number of bounces, is large. The relative im- 
portance of reflection and absorption losses can be deduced by looking 
at the ratio, x, of the quantities in the exponent in (72), i.e. 

x = _L cos 0,- | In R | . (73) 

ctsd 

For x > 1 reflection losses dominate, whereas for x < 1 the absorption 
losses are relatively more important. 

If we require a certain minimum transmissivity, (72) and (73) fix an 
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upper limit to the linear dimensions of the modulator in terms of a loss 
parameter. Thus if x > 1 and we set R lsd 00t ' = \ in (72), correspond- 
ing to a reflection loss of 50 per cent, we obtain 

- < * WV^ , n o (7 a) 

when 1 — R « 1. By comparing (71) and (74) we can establish the 
relative importance of diffraction effects and reflection losses in deter- 
mining the maximum allowed value of the ratio S/d* We find that the 
reflection losses are the important factor whenever the reflectivity is 
small enough that 

1-/2>0.17^-- (75) 

Am n 

and diffraction effects dominate when the inequality is reversed. For 
A = 1m, X« — 3 cm (X-band), the quantity 0.17 Xo/X m is of order 5 X 
10 -6 . Since even with the best of dielectric multilayers 12 it is difficult to 
exceed R = 0.999 (1 — R = 10~ ), it is evident that reflection losses 
are the important factor in materials for which k'/uq ^ 100, as is nor- 
mally the case. An important exception is a ferroelectric operated not 
far above its Curie point, for which the ratio k'/uq can be very large. 

If we now have the situation where x < 1 and we set exp ( — a s S 
esc 6i) = l /i in (72), corresponding to an absorption loss of 50 per 
cent, we obtain 

S<I-^l n2 (76) 

with no restriction on d. This equation then sets an upper limit on the 
modulator length when reflection losses are negligible. For materials 
having good optical quality where a s is typically of the order of 0.05 cm - f 
such a case might correspond to having R = 0.999 and d ~ 0.2 mm. 

5. 5 Effect of Lack of Parallelism in the Reflecting Surfaces 

A lack of parallelism in the mirrors causes the light and modulating 
waves to progressively slip out of synchronism, and so degrades the 
modulator performance. The effect may be estimated as follows. 

* By comparing (71) with (74) we are in effect determining the minimum allowed 
value of the ratio S/d(or a fixed sideband amplitude. Such a comparison is of sig- 
nificance since the sidebands carry the modulation information and since their 
amplitude will determine the signal-to-noise ratio at the detector. 

f Accurate values of the absorption coefficients of highly transparent materials 
are not readily available and are liable to depend on the details of individual sam- 
ple preparation. A value as ~ 0.05 cm -1 is a conservative estimate (cf. Ref. 13). 
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The time T taken for the optical signal to travel the length of the 
modulator is given by 



T 



Jo V«( 



dy 



v (y) 



(77) 



where v y = (c/n ) sin 6(y) is the y component of the optical phase 
velocity vector (Fig. 7). If the two reflecting surfaces are inclined at an 
angle xj/, the angle of incidence changes by \p with each successive bounce. 
When the total number of bounces N is large, we may approximate 6(y) 
by the continuous function 6(y) = 0, + (N^/S)y, so that the total 
optical transit time is given by 



'-7JM* + tF' 



dy 



Wo a • n i n o Nxf/S 

tt — S sin di -\ V- cos 6i , 

c c 2 



for 



Nf « 1 . (78) 



If the optical and microwave phase velocities are matched at the modu- 
lator input (sin 0,- = Wo/\/k')> * ne first term entering in (78) is just the 
microwave transit time. Thus at the output of the modulator the optical 
and microwave signals have shifted, in terms of the microwave phase, 
by an amount 

ntfam S 2 cos 2 0» 
c 2d sin 0i 

With the synchronism condition sin 0,- = n /vV, this becomes 

ttvVS 2 



A<p = 



d\. 



(1 - no 2 A')*, 



(79) 




/ g(0) 1 
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Fig. 7 — Model used in calculating the effect of nonparallelism of the reflecting 
surfaces. 
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where X m is the free-space modulating wavelength. This phase shift 
must be less than t if we are to avoid stripping modulation off the beam, 
since the effect of the modulating field will reverse as the phase shift 
exceeds r. This leads to a condition on the parallelism of the mirrors 

d\ m 1 , _* 

■Vk'S (1 — Wo A ) 

5.6 Effect of Strains 

A complete evaluation of the effect of strain on the modulator per- 
formance is very difficult; hence we shall not attempt a complete treat- 
ment here. Nonetheless, since strains can have a profound effect on the 
operation of the modulator, even an approximate treatment such as 
presented herein can be of value. 

In the present discussion it is convenient to consider two broad classes 
of strain, uniform strain and random strain. The effect of a strain uni- 
form throughout the crystal can be described in terms of a photoelastic 
matrix 14 which must then be combined with the electro-optic matrix in 
order to completely describe the optical properties of the crystal. The 
net effect is to add a static birefringence to that induced electrically. 
This problem can in principle be solved analytically, although a large 
strain produces an additional complication since it will in general lower 
the crystal symmetry and hence change the form of the electro-optic 
tensor. Although cases of quasiuniform strain sometimes occur,* we shall 
consider here only the case of random strain, since this imposes a more 
fundamental limitation on device performance. 

In a randomly strained crystal where the strain-induced variations 
occur over distances large compared with an optical wavelength, a ray 
in progressing through the crystal will see a continuously varying index 
ellipsoid. If there is little correlation between the shape and orientation 
of the index ellipsoid at two adjacent points, then a statistical model may 
be used to study the effect of the strains. Such a calculation is now being 
carried out and will be reported separately. We confine ourselves here to 
some general remarks, and a simplified analysis which is valid in the 
limit of very severe strain. 

The general effect of random strain is to partially destroy the spatial 
coherence of a plane wave traversing the medium, without affecting its 



* For example, cubic crystals drawn from a melt are sometimes observed to be 
uniaxial, with the growth axis as the optic axis. 
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temporal coherence.* The wave phase at the output of the modulator 
may then be written as the sum of a coherent and incoherent term, i.e. 

*(r,0 = <p{t) + S(r), 

where 5(r) is a random variable, and r is some suitably defined transverse 
position vector. Random strains will also mix the two normal modes of 
the unstrained crystal. We shall not consider the details of this mixing 
here, but simply assume that there exist at the output of the modulator 
two orthogonally polarized modes whose wave phases are given by 

*i(r,0 = <pM + «x(r) (81) 

and 

*ii(r,0 =*i(0 + *i(r) (82) 

where <p x and <p% are the wave phases in the unstrained crystal as previ- 
ously calculated, and 8 ± (r) and 5 B (r) vary more or less randomly across 
the surface of the beam. In this model the effect of strain on the modula- 
tor performance depends on the autocorrelation functions of 8 X and d D , 
and the cross-correlation between 5 ± and 8 t , as the subsequent analysis 
will make clear. These correlations will depend on the state of strain of 
the crystal, and on the length of the optical path, becoming less as the 
path length increases. 

The electric vectors at the output of the modulator can be written in 
the form 



e ± (i,t) = E ± exp {j[wt + Vx (t)] + j8 x (r) } (83) 



and 



e,(r,0 = E t exp [j\wt + „,(*)] + j«,(r)} (84) 

where E ± and E\ are the wave amplitudes. From these equations we can 
show that although strains do not affect the basic phase modulation 
process, they do affect the over-all modulator-detector system perform- 
ance. Details depend on the system used, but the following general state- 
ment can be made: Any optical communication system in which the 
modulation or detection process is based upon the heterodyning or homo- 
dyning of two optical signals requires that the two signals have relative 
spatial coherence. Any process that tends to destroy that coherence, such 
as random strain, will degrade the system performance markedly. It is 

* The fact that in general strain will also degrade the collimation of the beam 
is of importance in some applications. See Appendix B.3. 
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therefore necessary to be able to evaluate which system will allow us to 
recover the modulating intelligence with maximum certainty. All modu- 
lation-detection schemes presently under consideration fall into three 
classes 

(i) amplitude modulators followed by quantum counters, 
(ii) phase or amplitude modulators followed by optical heterodyne or 
homodyne detectors, and 

(Hi) phase modulators followed by optical frequency discriminators 
and then by quantum counters. 

Starting with (83) and (84), we analyze these three classes in Appendix 
B. The results are summarized below. 

5.6.1 Amplitude Modulators followed by Quantum Counters 

As discussed in Section IV, the amplitude modulator consists of a sys- 
tem in which two phase-modulated waves of the form given by (83) and 
(84) are combined in an analyzer. As shown in Appendix B, the time- 
varying intensity transmitted by the analyzer is proportional to an in- 
tegral of the form [see (118)] 



L 



sin 

B COS 



5 ± (r) - «,(r)ldo (85) 



the integration being taken over the cross section of the beam. In the 
case of severe random strain, when 8 X and 8] can independently take on 
all values (i.e., no cross-correlation between 5^ and S\\ ), the integral (85) 
vanishes and there is no amplitude modulation. More moderate strain 
would result in a nonzero depth of modulation which would, however, 
be less than that attainable in a strain-free crystal. 

5.6.2 Heterodyne or Homodyne Operation 

In heterodyne or homodyne operation the modulated signal is mixed 
with a second locally generated optical signal to obtain a microwave in- 
termediate frequency which is then detected conventionally. Assuming a 
square-law detector and a strong local oscillator having both spatial and 
temporal coherence, the output signal at the difference frequency is pro- 
portional to [Appendix B, (121)] 



f cos 

Js sin 



|*x.i(r)|do (86) 



which vanishes in the limit of very bad strain where 8 X and fig can take 
on all values. 
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Thus in this mode of operation, too, strains result in a reduced detector 
output. The effect here is likely to be more severe than in the case of the 
simple amplitude modulator-quantum counter system, since for a given 
degree of strain the autocorrelation of 5 x or S B is expected to be less than 
the cross-correlation of 8 ± and 5 D . Thus, in a given crystal the integral 
(85) is liable to have a larger value than the integral (86). 

5.6.3 Phase Modulators Followed by Frequency Discriminators 

In the third system under consideration phase modulation is directly 
converted to amplitude modulation in some frequency-selective optical 
circuit 15 ' 16 of which a Fabry-Perot etalon is perhaps the simplest example. 
Spatial coherence per se plays no essential role here. Analytically, as is 
shown in Appendix B, the integral analogous to the ones given by (85) 
and (86) which arises in this case has the form [see 1241 

flexp [j8 ± ,,(r)]\ 2 da. (87) 

Js 

This of course has a constant value equal to the beam area and can 
therefore not influence the system response. However, strain can still 
adversely effect this system to the extent that the discriminator per- 
formance is degraded by loss of collimation in the transmitted beam. 

Thus strain will always be a problem in optical modulators. In fact the 
preceding considerations show that the presence of severe random strain 
will result in no useful output from those modulation-detection schemes 
which rely on the interference of two beams. In this case only the phase 
modulator-frequency discriminator system has any real chance of suc- 
cess. 

VI. OPTIMUM DESIGN PROCEDURE 

6.t Choice of Modulator Dimensions 

For the purpose of discussing the optimization of modulator dimen- 
sions as well as the choice of operating temperature, a useful figure of 
merit is the product of the modulation sensitivity and the optical trans- 
missivity, which for all modulators can be written as 

* i _ ^ „ .« . rsi fci . (•& , (** 



8* 



VF 
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Here d is the distance between the reflecting surfaces (either h or w, de- 
pending on whether the 0° mode or the 90° mode is used), and I is the 
width of the reflecting surfaces; recall that S is the modulator length. 
The remaining functions f(a m S/2) and p(n /\/7') have been defined 
above in (12) and Tables I and II respectively. Equation (88) applies 
when we have velocity matching. 

The optimum modulator dimensions can now be found by maximizing 
the figure of merit. It is evident that the dimension I should be made as 
small as possible. If the beam diameter D b is specified, this then deter- 
mines the minimum value of I, l m - ia = D b . The mirror spacing d may 
next be chosen so as to maximize (88), the appropriate value being 

d opt = -25 In R ?-<- (l - *£Y. (89) 

The ratio S/d obtained from this equation corresponds to a reflection 
loss of roughly 50 per cent. It is worth pointing out that for practically 
attainable reflectivities, this optimum value of S/d is less than the maxi- 
mum value set by diffraction effects [see (77)], and so represents a reason- 
able design value. 

When the optimum value of d given by (89) is substituted back into 
(88), there results an equation for g which depends only on S; S may 
then be chosen so as to maximize this figure of merit. The optimum 
value, S op t , is defined by the transcendental equation 

1 1 \Q a a 



, > ■ (90) 

exp ttmOopt — 1 AOlmOopt ^0 «m 

In the absence of optical absorption (as = 0), the optimum length is 
given by a m S opt = 1.25. As a 8 increases, the optimum length decreases, 
as shown in Fig. 8. When the inequality (n /\/ r K , )(a m /a B ) « 1 is satis- 
fied, we find that the optimum length is determined solely by the optical 
attenuation, and is given by 

The dimensions of the modulator are thus completely determined: I 
by the beamwidth, d by the optical reflection losses and S by the micro- 
wave and optical attenuation constants. For materials having sufficiently 
low microwave and optical losses, the optimum length given by (90) 
may be greater than any practically realizable crystal length. Under 
these circumstances one should make S as large as possible, and then 
choose d to satisfy (89). 
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Fig. 8 — Dependence of optimum modulator length on optical and microwave 
loss parameters. 

6.2 Choice of Operating Temperature 

The dielectric and electro-optic constants of useful modulating ma- 
terials are strongly temperature-dependent. In ferroelectrics, for exam- 
ple, both the electro-optic coefficients and the microwave losses increase 
dramatically near the Curie point. We should thus like to establish 
whether there exists a preferred operating temperature which maximizes 
the figure of merit [see (88)] independently of the modulator dimensions. 

We begin by assuming that the modulator dimensions have been op- 
timally chosen according to the results of the preceding section. With 
this choice of dimensions, (88) becomes 



-> 2tt 3 

^opt = — n 
Ao 



V?D b In R 



h r *i 



l* 



/ no \ r no/V? T 
°eJ P W2) Ud - *PA0'J (92) 



X e 



-<[v^/g 



where u = a m S opt is found from (90) and where we introduce the nota- 
tion ( a a /a m )(V« 7 /no) = Q. The function g(u) = \/ue~ Qu f(u/2) is, 
through (90), an implicit function of the independent variable Q. The 
function may be approximated as 
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[l-exp^W)] 
g(Q) » 0.214 U -^ LA . 

For a TEM-like wave, the microwave attenuation constant is given in 
terms of the dielectric constant and loss tangent of the medium by 

a m = ~ Vk' tan 5. (93) 

Am 

The figure of merit can now be approximated solely in terms of material 
parameters as 

nn 4 i 1 IW) 

gopt = 15.8- ^/_ asDblnRK , 

[pE0> (94) 



AsM* 



,1 - exp| -- 



f _ 3 / «sX OT y 
|_ 2 \7rn0 tan 6/ _ 



(1 - nf/W ~ J 

where all the numerical constants, including the impedance of free space, 
t) = 377 ohms, have been combined into a single factor. 

Once the temperature dependence of the electro-optic coefficient, di- 
electric constant and loss tangent are known,* (94) can be used, at least 
in principle, to select an optimum operating temperature. 

Further analytical progress requires some assumption of the tempera- 
ture dependence of these coefficients as well as a choice of the mode of 
operation. To be definite we select the mode for which p(rio/W) = 
s/k'/uo . As may be seen from Tables I and II, this is the most favorable 
value that can be realized and applies, for example, to the linear, 0° phase 
modulator. 

The temperature coefficient of the linear electro-optic coefficient is 
expected to follow that of the dielectric constant. 1 The reason for this is 
as follows: The electro-optically induced birefringence arises funda- 
mentally from a change in the polarization, rather than the applied field, 
and in fact Pockel's 18 original formulation of the electro-optic effect was 
in terms of constants e,y relating the coefficients of the index ellipsoid to 
the induced polarization. The r,y coefficients used here are related to 
Pockel's etj by 

Tij = € (k — !)««/. 



* The index of refraction, no , to the extent to which it is determined by the 
electronic polarizability, should be almost independent of temperature. 
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Thus for materials having a high dielectric constant we expect r^ to be 
approximately proportional to k', and to the extent that the coefficients 
e,, are temperature independent, 17 to have the same temperature depend- 
ence as k. 

A similar argument in the case of the quadratic effect leads to the con- 
clusion that p should be proportional to k' 2 . Here, however, the situation 
is complicated by the fact that the bias field E , which for optimum op- 
eration should be as large as the breakdown field strength of the material 
permits, will be temperature-dependent whenever this breakdown field 
strength is temperature-dependent. In many solids we expect dielectric 
breakdown to occur at some definite value of internal field. In materials 
of high dielectric constant the internal field is determined by the polari- 
zation, so that the breakdown field strength, and hence the allowed dc 
bias, will be inversely proportional to «'. The net result is that the product 
of pE may be taken proportional to the first power of k', and the quad- 
ratic and linear effects become formally identical. 

The preceding assumptions allow us to write the figure of merit for a 
given material wholly in terms of its dc and optical dielectric parameters, 
so that $ becomes proportional to 



k"' 



(1 - no 2 A') 1 



D-- »(-s££j]- (95) 

We must now distinguish between low and high microwave losses, the 
criterion being the magnitude of the ratio 3\ m a s /2irn tan 8. For suffi- 
ciently low microwave losses, the exponential factor in (95) will dis- 
appear, so that the figure of merit is proportional simply to k' 3/2 / 
(1 _ no 2 A') 1/4 - Since the ratio n 2 A' will always be substantially less 
than unity, the optimum operating temperature is that which maximizes 
k'. For ferroelectric materials this will mean lowering the temperature to 
the vicinity of the Curie point, T c . Of course the losses will normally 
increase as the temperature is lowered, so that at some point the approxi- 
mation of low microwave loss breaks down. For sufficiently high loss 
we may expand the exponential term in (95), obtaining a figure of merit 
proportional to 

1 A /JE_ ( 96 ) 

(1 - n 2 A') i \ tan 5' 

For ferroelectrics such as KDP 19,20 or SrTi0 3 21,22 in which both k' and 
tan 5 are dominated by a (T — T c )~ l dependence, T c being the Curie 
temperature, this result says that operation near the Curie point is 
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favored (although not strongly) even in the presence of substantial 
microwave loss. 

6.3 Bandwidth Limitations: Design of Optical Coupler 

The question of the bandwidth of the zigzag modulator has two as- 
pects — microwave bandwidth and optical bandwidth. Since both the 
optical and modulating signals propagate as TEM-like waves, it might 
at first be thought that the bandwidth would be limited only by the dis- 
persion in the material parameters of the modulating medium. For most 
materials this would imply a microwave bandwidth as large as 10 or 20 
gc, and since the fractional optical bandwidth required to accommodate 
even microwave frequency sidebands is very small (of order 10~ ), op- 
tical dispersion should never be a problem. 

A more discriminating analysis of the question of bandwidth reveals, 
however, that the physical dimensions of the modulator limit the realiza- 
ble bandwidth, especially when these dimensions are optimized so as to 
obtain the ultimate in low-power performance. In an optimum design 
the length <S opt and width (or height) d op t are uniquely determined by the 
electrical properties of the modulating medium both at microwave and 
optical frequencies as well as the loss properties of the medium [see (89) 
through (91) and the attendant discussion]. For highly transparent low- 
loss microwave materials, d opt computed from (89) is nominally in the 
order of a few tenths of a millimeter and S opt is usually several centime- 
ters (see Table IV in Section VII for typical results). The one dimension 
which remains unspecified is the beam diameter Db . As remarked above 
in Section 6.1, Db should be chosen as small as possible and yet remain 
within the diffraction limits of the device over the length S opt . Ac- 
cordingly, we find that, based upon certain operational practicalities 
and diffraction limits, Db must normally be restricted to about a milli- 
meter. Thus, as the aperture of the device is typically an order of magni- 
tude larger than the spacing between mirrors, it is impossible to satisfy 
(70); i.e., the beam bounces back out of the modulator, and hence can- 
not be confined to the active area of the device by means of a simple 
aperture arrangement. In such cases the beam can be confined to the 
interior of the modulator only if we include an appropriate optical coupler 
at the input (and output) of the modulator. 

It is shown in Appendix C that a suitable coupler which allows all of 
the incident optical energy to be transferred into the active medium of 
the modulator is a Fabry- Perot resonator of the form shown either in 
Fig. 10(a) or Fig. 11. The physical operation of the coupler is easily un- 
derstood when one observes that because of its finite length energy must 
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leak off to the right into the interior of the modulator. Then by properly 
adjusting the spacing between reflectors as well as the reflectance of the 
lower reflector it is possible to completely cancel the reflected wave, and 
have the energy trapped in the coupler section walk off into the modula- 
tor. Expressions for the required spacing and reflectance are derived in 
Appendix C. Here we simply point out that such a coupler can be realized 
in a typical case with reflectivities in the neighborhood of 99 per cent, 
which is easily attained, and with dimensions L c and l c of the same order 
as D b and d opt , respectively. 

We are now in a position to examine the bandwidth limitations im- 
posed by the coupler and modulator dimensions. The bandwidth is gen- 
erally determined by three factors. In order of importance these are 
(*') the spacing l c and loaded Q of the coupler, 

(it) the size of D b compared with the microwave modulating wave- 
length, and 

(to) the microwave impedance of the structure. 
The last factor, microwave impedance, is of practical significance be- 
cause the impedance of the modulator will inevitably be low, and it be- 
comes progressively more difficult to obtain a broadband match into a 
transmission line as its impedance is lowered. However, a number of 
standard microwave techniques are available to overcome this difficulty. 

The second factor, the beam diameter, affects the bandwidth in the 
following way. In order that the microwave modulation be detectable 
in a direct simple way it is necessary that the modulation frequency phase 
be approximately constant across a given section of the beam. Hence it 
is necessary that the projection of the beam cross section along the axis 
of the structure be appreciably less than half the microwave wavelength 
in the medium. Mathematically 

Db«o T— 7=7 (97) 

2 cos 9i Vk 

where \ M is the free-space wavelength. 

Finally, with reference to the first factor, the coupler spacing and 
loaded Q, we note that in order to avoid microwave impedance discon- 
tinuities in the low-impedance modulator line it is mandatory that l c and 
d be very nearly equal. If the spacing d is then kept small so as to obtain 
the lowest possible modulating power, the input coupler will limit the 
bandwidth in the following way. The Fabry-Perot coupler will be res- 
onant at a series of frequencies spaced by 

A/= / - . (98) 

2n l e cos 6i 
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The base bandwidth 8f of the coupler will be determined by the Q of an 
individual mode, and can be shown from Appendix C to be 

. „ fVl - n 2 sin 2 di for n 2 < \/k' 

IT Db no \ n c 2 . /—, 

(cos di for n > vr< 

Since this bandwidth is available in each of the Fabry-Perot modes, 
microwave subcarrier operation at frequencies above this baseband is 
allowed 23 provided that the subcarrier frequency (or optical sideband) is 
some multiple of the mode spacings given by (98). The problem, of 
course, is that the mode spacings become very large when the coupler 
reflector spacing l c is small, and since l c must be comparable with the 
modulator thickness d, it follows that the lowest allowed microwave sub- 
carrier will be at much too high a frequency when d is adjusted to its 
optimum value of a fraction of a millimeter (e.g., d = 0.2 mm cor- 
responds to A/ ~ 300 gc). It is evident then that the minimum power 
configuration is restricted to baseband operation with a bandwidth de- 
termined either by (97) or (99), whichever is smaller. For Db = 1 mm, 
(97) gives f m < 10 gc when the modulating medium is a ferroelectric 
such as strontium titanate, while (99) gives 8f tt 5 gc. 

Greater bandwidth can be obtained, at the expense of increased modu- 
lation power, by making the thickness d considerably greater than opti- 
mum. If d is made sufficiently large (e.g., 3.6 mm in the case of the stron- 
tium titanate modulator), a simple input aperture can be used in place 
of the Fabry-Perot coupler.* This removes the frequency limitations set 
by (98) and (99), leaving only that imposed by the finite beam width 
(97) and material dispersion. The frequency limitation imposed by the 
beam width is quite fundamental. In fact, when the projection of the 
beam diameter along the axis of the modulator is exactly equal to a full 
microwave wavelength in the material there results no net modulation. 
We therefore define a cutoff frequency for the modulator as 

'• - fto-av? - (100) 

In summary, then, we find that in order to fully exploit the low modu- 
lation power requirements of the zigzag modulator it is necessary to re- 
strict the frequency of operation for most substances to a few gc at base- 
band,! and to use a relatively sophisticated coupler to get the optical 

* If no 2 > Vk' it will be necessary to use a more complex aperture analogous 
to the one shown in Fig. 11. 

t For CuCl the operating frequency can be as high as about 20 gc (see Table IV 
in Section VII). 
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energy into the device. By relaxing the power requirements the upper 
frequency limit can be extended above 10 gc in a relatively simple struc- 
ture. Examples of both kinds of design are given in the following section. 

VII. DISCUSSION AND CONCLUSIONS 

We have analyzed in detail the behavior of zigzag phase and amplitude 
light modulators (Sections III and IV) as well as the basic underlying 
theory for this class of modulators (Section II). Our treatment differs 
from the work of other investigators'- 3 along these lines in a number of 
crucial ways. We have included microwave losses in the initial calculation 
of the optical retardation and have considered the 90° mode of operation 
in addition to the originally proposed 0° mode for several different classes 
of electro-optic media. The analysis has been detailed enough and has 
treated enough of the important practical factors which limit the opera- 
tion of the device (Section V), so that a realistic optimum design pro- 
cedure can be presented (Section VI). In essence, this procedure requires 
that the dimensions of the modulator be so chosen that all of the incident 
microwave power be confined to as small a volume of the electro-optic 
material as the microwave, optical, and reflection losses will permit. As 
a result, large modulating field strengths are obtained with moderate 
amounts of power. 

Based upon our treatment we find that there are three main aspects 
to broadband modulator design. The first is the velocity synchronism 
condition; here we find that the basic phase velocity matching criterion 
is unaffected by microwave loss. The second is the effect of modulator 
dimensions and aperture on bandwidth. These problems have been dis- 
cussed in Section 6.3, where it was shown that the aperture limits the 
highest modulating frequency in a fundamental way [see (100)] and that 
the other dimensions in effect determine the requisite coupling scheme 
(either Fabry-Perot or simple aperture) and thence the bandwidth. The 
final criterion of significance is the strain properties of the modulator 
crystal (Section 5.6). Crystalline strain is a serious problem in all forms 
of electro-optic amplitude modulators, but is particularly deleterious in 
the traveling-wave modulator because of the long optical path. 

We. are now in a position to present some actual designs of zigzag phase 
and amplitude light modulators as well as the performance to be ex- 
pected of these devices. Before doing so, however, we first summarize 
in Table III the pertinent microwave and optical material parameters, to 
the extent that they are known, for a number of electro-optic modulating 
crystals which are currently available. With regard to the data presented 
in Table III, much of which is tentative, two points should be noted. 
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The first is that the microwave properties for the different materials 
listed have not all been measured at the same frequency. For simplicity 
we shall assume that there is negligible dispersion in the microwave 
parameters. The same assumption is made for the optical data. The sec- 
ond point to be noted is that the data available for the electro-optic co- 
efficients apply in most cases at very low frequencies. Hence the values 
for the electro-optic coefficients represent the undamped values and 
consequently overestimate in many cases the desired microwave or 
clamped values for these coefficients. Thus this table should be viewed 
simply as a guide to the material parameters. From the data given in 
Table III we select three materials, CuCl, SrTi0 3 and KDP, as being 
the best in their respective crystalline classes. Each of these materials 
combines good microwave properties (low loss with high intrinsic im- 
pedance) with a large electro-optic coefficient. With reference to SrTi03 
it should be noted that the quadratic effect exhibited by this material is 
anomalously high, as it is with BaTi0 3 and KTa0 3 , the effect being 
10 5 -10 7 times higher than in the best glasses. When one couples this to- 
gether with the fact that the "effective linear coefficient" for quadratic 
materials is pE , one immediately sees from the data listed in Table III 
that for E > 1000 v/cm, quadratic materials are potentially as good as 
if not better than linear materials. 

It is instructive to carry out two designs. One, based on the optimum 
design procedures of Section 6.1, results in a modulator requiring an 
absolute minimum of modulating power. However this design requires 
the use of an input optical coupler and is inherently narrow-band. In 
another, more practical design, the mirror spacing d is chosen to allow 
the use of a simple aperture. In this design we obtain increased band- 
width at the expense of higher modulating power. Tables IV and V sum- 
marize the design parameters, expected performance, and several of the 
practical limitations for each of these designs. Note that we do not in- 
clude an amplitude modulator design for KDP. The reason, of course, 
is the difficulty encountered in trying to compensate for the static re- 
tardation in such a material (See Section 3.3.) 

In order to carry out the calculations several design parameters were 
chosen independently as a matter of convenience. Thus the optical wave- 
length was taken as the visible red fine of the He-Ne gas laser at 6328 A, 
and in the absence of experimental data the optical absorption coefficient 
a a was assumed to have the same value as that found in good optical 
glasses, viz., a s ~ 0.05 cm -1 . The multilayer reflectivity was taken as 
0.999, about as high as is practically achievable. 12 A reflectivity this high 
is essential if the ultimate in low-power modulation is to be obtained in 



2662 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1963 

the optimum design, but the requirements can be relaxed considerably 
in the practical design. For example, a decrease in reflectivity to 0.982 
increases the optical insertion loss of the practical SrTi0 3 modulator by 
only 1.7 db. 

In the absence of any really complete temperature data we have taken 
room temperature, 300°K, as an operating temperature. Although it is 
true that for the ferroelectrics, SrTi0 3 in particular, the temperature 
analysis of Section VI shows that the modulating power is reduced some- 
what as we approach the Curie point, the concomitant increase in k 
reduces 0; below the practical minimum and at the same time reduces 
the available bandwidth by decreasing the microwave wavelength in the 
material [see (97) and accompanying discussion]. Thus 300°K is in fact 
close to the practical optimum operating temperature for SrTi0 3 . 

The operating frequencies given in Tables IV and V have been arrived 
at by the following considerations. A practically realizable aperture 
which remains within diffraction limits is nominally 1 mm. The optimum 
design therefore requires the use of a coupler with restricted bandwidth. 
The data of Table IV are thus calculated for a modulating frequency of 
5 gc, which is within the passband of all modulators. In Table V the 
operating frequency has been chosen as 10 gc and is in all cases less than 
the cutoff frequency. What should be noted here is that by ignoring the 
bandwidth problem, we obtain in the optimum design shown in Table IV 
the irreducible minimum microwave modulating power. In the case of 
the practical design listed in Table V the spacing between mirrors d has 
been increased to allow the use of a simple aperture. The price we pay, 
of course, is in increased modulating power. We should also make special 
note of the fact that from a practical viewpoint the impedance level of 
the design given in Table V is much more favorable than for the one 
given in Table IV. The restriction on the maximum attainable modulat- 
ing frequency in the practical design is the aperture size [see (97)]. 

The column headed "configuration" in Tables IV and V deserves some 
comment. For each material considered, the configuration chosen was 
the one for which the function p(ft /vV) had its maximum value. In 
the case of SrTi0 3 there are three possible configurations for the phase 
modulator. Table I shows that the 0° phase modulator (PM) in the 
_L-mode and the 90° PM in either the _L or ||-modes are all equally 
suitable choices. We should bear in mind at this point that for SrTi0 3 
(as well as BaTi0 3 and KTa0 3 ), it is only the difference in the electro- 
optic coefficients | pi — p 2 | that is known (see Table III) and not p t 
or p2 alone. However, since in an amorphous material p 2 is identically 
zero (Appendix A) we expect that in a cubic material | p 2 | will still be 
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substantially less than | pi |, so that we may take the measured value of 
I Pi — Pi | as an approximate measure of | pi | alone. This has been done 
wherever necessary. 

The remaining parameters listed in Tables IV and V require little or 
no further comment. The angle of incidence 0,- necessary to achieve 
velocity synchronism is given together with the deviation in 0, , | A0, |, 
which results in a 3-db degradation in modulation sensitivity. The angle 
^max is the maximum permissible angular deviation in the parallelicity 
of the mirrors if the performance is to be degraded by less than 3 db. 
The modulator dimensions are contained in Db , d opt or d, and »S op t or S 
depending upon which design is being considered. The column headed 
IL in the tables refers to the optical insertion loss suffered in traversing 
the modulator. Finally, the quantity P(M = \) is the microwave modu- 
lating power required to achieve a modulation index of 0.5. 

It is worth mentioning at this point that Tables IV and V implicitly 
contain several idealizations. The first is that the crystals are strain-free. 
For KDP this assumption is approximately justified, but for CuCl and 
SrTi0 3 it is far from reality and, as noted previously, could seriously de- 
grade the performance of the amplitude modulator. The second idealiza- 
tion comes from the assumption that the modulating field is uniform 
over the entire cross section of the crystal. In any practical modulator 
this will not be strictly true. However, to the extent that fringing fields 
can be neglected in a parallel plate line completely filled with electro- 
optic material, this assumption is realizable (see Section 5.1). In conse- 
quence, the calculations made for the required modulating power should 
be reasonably accurate. The primary factor which would cause an in- 
crease in modulating power is a deviation from the correct angle of in- 
cidence. However, the tolerances in 0, given in Tables IV and V are 
sufficiently large that maintaining the required angle should present no 
problem. Finally we should point out that we have assumed that a dc 
bias as large as 30 kv/cm can be applied to the SrTi0 3 sample without 
dielectric or interfacial breakdown, and that the electro-optic effect does 
not saturate at these high fields. It is clear that a high dc bias field is 
desirable with the SrTiOj modulator, since this increases the effective 
linear electro-optic coefficient pE . 

The results given in Tables IV and V demonstrate that broadband 
microwave modulation of light can be effected with very little modulat- 
ing power. In a practically realizable design, we find that a 1 mm X 3.6 
mm X 1.4 cm rectangular slab of SrTi0 3 biased with 30 kv/cm should 
produce 50 per cent linear modulation of light with less than 5 watts of 
power at any modulating frequency from dc to 10 gc. A CuCl slab of 
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dimensions 1 mm X 0.75 mm X 6.7 cm should require still less power 
(«l-2 w) without the need of any dc bias. At the expense of reduced 
bandwidth and increased optical complexity, this power can be reduced 
to less than a watt by a correct choice of the modulator dimensions (op- 
timum deisgn). In the zigzag configuration, CuCl appears to be slightly 
better than SrTiOa , in the sense of achieving moderately broadband 
performance with low modulating power. Cuprous chloride has the ad- 
vantage of extreme bandwidth, of not requiring any bias, and of not 
needing reflecting mirrors since 0,- exceeds the critical angle for total re- 
flection. Compared with cavity-type modulators using KDP or ADP 
the CuCl zigzag modulator typically requires three orders of magnitude 
less drive power and has up to two orders of magnitude greater band- 
width. Compared with other traveling- wave systems, 3 or with those 
iterative structures proposed to date, 5 it is substantially superior from 
the viewpoints of size and power. 

As far as the practical realization of these zigzag modulators is con- 
cerned, we believe that there are no insurmountable constructional prob- 
lems. However, it should be noted that at the present time single crystals 
of suitable cubic electro-optic materials of the required optical quality 
and size are not readily available. Continuing materials research, how- 
ever, indicates that this situation should improve. 

We should point out that it appears that high dielectric constant cen- 
trosymmetric paraelec tries are particularly well suited to the zigzag con- 
figuration. This stems from the fact that the dielectric constant of these 
materials is intermediate between those of conventional dielectrics and 
those of ferroelectrics, and their quadratic electro-optic effect is large. 
The dielectric constant is large enough in these materials that the angle 
required by the synchronism condition is reduced to the smallest prac- 
tical value. This results in a long optical path in a crystal of modest 
physical size, and hence in very efficient use of modulatng power. The 
existence of a large quadratic electro-optic effect allows us to obtain a 
very large effective linear coefficient by using a suitable dc bias. It is of 
interest to note in this connection that the performance of the SrTi0 3 
modulator can be improved somewhat by going to lower temperatures. 

In closing, we observe that the volume of active material required in 
the zigzag modulator can be cut in half by using the double-pass ar- 
rangement shown in Fig. 9. Here both the light and the modulating sig- 
nal are reflected from the end of the modulator, with velocity synchron- 
ism being maintained for both incident and reflected waves. In this way 
the physical length of the crystal required to achieve a given electrical 
length can be halved. 
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Fig. 9 — Double-pass zigzag modulator configuration. 
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APPENDIX A 



In this appendix we derive expressions for the electric field dependence 
of the refractive indices of the two normal modes used in the zigzag 

modulator. 



A.l Linear Electro-Optic Effect in Cubic Materials of Symmetry Ta 

For this class of materials, to winch cuprous chloride belongs, the equa- 
tion of the index ellipsoid relative to the crystalline axes x',y',z' is 

— 2 Or' 2 + y n + z' 2 ) + 2r ii (E z ,y'z' + E v ,z'x' + E t ,x'y') = 1 (101) 

where n is the index of refraction in the absence of the electric fields 
E x i , E y i , and E z i , and r« is the electro-optic coefficient. For an electric 
field applied along the z' axis (0° modulator) in a crystallographic [001] 
direction, the equation reduces to 



—Ax + y + z ) + 2niE,ixy = 1 
n 2 



(102) 



which is the equation of an ellipsoid having one of its principal axes along 
the z' axis, and the remaining two in the x'y' plane at 4.5° to the x' and 
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y' axes.* In a coordinate system x,y,z coincident with these axes the 
index ellipsoid has the normal form 

(L + Ui e)x 2 + {\- r a E)y* + £ = 1. (103) 

\w 2 / W / n 2 

The principal indices of refraction and their respective crystallographic 
directions are thus 



along [110] 

ir,,,sr «.; - , — ^V ,, along [TlO] 
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(104) 



In the 90° modulator, the field is applied in the .r-direction. For this 
case one principal axis of the index ellipsoid lies along x', and the re- 
maining two in the y'z' plane, at 45° to y' and z'. Again, taking these 
principal axes as coordinate axes, we find for the principal indices of 
refraction 



90° case 



2 2 

n x = no 


along [100] 


2 Wo 


along [011] 


Uy 1 - Wo 2 r 4 i#* 


2 Wo 

n, = ■—, = — ^r 


along [Oil] . 



(105) 



1 

A.2 The Quadratic Electro-Optic Effect in Cubic and /or Amorphous Mate- 
rials 
All substances (solids, liquids, and gases) which have a center of 
symmetry exhibit the quadratic Kerr effect. For most amorphous sub- 
stances the electro-optic Kerr effect may be described by the empirical 
relation 

An = K\ Q E\ (106) 

In this relation An is the induced birefringence, K is the Kerr constant, 



* It is not essential that the electric field be applied in a [001] direction. In the 
general case, however, a principal axis of the index ellipsoid will not be in the di- 
rection of the applied field, as we assumed in our analysis in Section II. This com- 
plicates both the analysis and operation of the modulator, and results in a 
decreased electro-optic effect. The basic approach used here would nonetheless re- 
main valid, since for cubic materials the orientation of the index ellipsoid depends 
only on the direction of the applied field, and not on its magnitude. Thus, there is 
no "wobbling" of the index ellipsoid, and the two normal modes are not mixed 
in the modulation process. 
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X is the free-space optical wavelength, and E is the applied field. Thus 
in the presence of an electric field isotropic materials effectively become 
uniaxial and will therefore display the phenomenon of double refraction. 
For crystalline solids the effect is more complicated, and must in 
general be described in terms of a matrix of coefficients. Thus in the 
general equation for the index ellipsoid 

a n x' 2 + a^f + a 33 z' 2 + 2a n x'y' + 2a 23 y'z' + 2a 31 z'x' = 1 

the coefficients have the field dependence 

flu - — = pnE x " + piiEy' + p\ 3 E z - + puE X 'E u > + puEy>E,> 

n* 2 (107) 

+ pnEfE* etc. 

leading to a 6 X 6 matrix of 36 independent electro-optic coefficients 
Pij . As in the linear case, however, crystal symmetry limits the number 
of nonzero coefficients, and in fact for cubic materials having the highest 
symmetry (Oh) the quadratic electro-optic effect matrix contains only 
three independent coefficients and is of the form 





Pij = 



Pi 


pi 


pi 
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(108) 



Note that for amorphous materials p 2 = 0, and p } = p 3 . 

If the field is now applied along a cube edge, e.g., the z' axis of a crys- 
tal of symmetry Ok , the equation for the index ellipsoid reduces to 



Too 2 



(x' 2 + if + z n ) + EJUx 2 + y n ) + Pl z' 2 ] = 1 . (109) 



This equation is in normal form as it stands and shows that the field 
induces a change in the refractive index isotropically in the x'y' plane 
as well as along the field direction or z' axis. The principal indices of re- 
fraction determined from (109) are 



0° case 



= n« = 



Too 
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For a field applied in the x direction (90° modulator) the correspond- 
ing results are 

n x 2 = - , n \ = , along [100] 

90° case (HI) 

along [010], [001]. 



2 2 no 

n„ = n, = 



1 + n<? P2 E x 2 

A. 3 Linear Electro-Optic Effect in Uniaxial Materials of Symmetry D u 

For an electric field applied along the c axis of a uniaxial crystal (taken 
to be the z' axis, so that we have the 0° case), the equation for the index 
ellipsoid is 

- 2 ix n + y") + - 2 z' 2 + 2*A*Y = 1 (112) 

where n and ?i e are respectively the ordinary and extraordinary re- 
fractive indices, and r 63 is the electro-optic coefficient of interest. As was 
the case for cubic materials, one of the principal axes of this ellipsoid 
lies in the crystallographic [001] direction, the other two in the [110] 
and [110], so that in a coordinate system x,y,z coincident with these 
axes the principal indices of refraction are 



2 no 
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0° case n* -- ^_tt along [TlO] (113) 

1 - no-mEz 
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For the 90° case the modulating field is still applied along the c axis, 
but this is taken in the x direction so that the principal indices of re- 
fraction are 
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APPENDIX B 



In tli is appendix we develop some expressions for the output of the 
three modulator-detector systems enumerated in Section 5.6, using (83) 
and (84) to account for the effects of random strain. These equations for 
the wave amplitudes at the output of the modulator were 

e ± (T,t) = E x exp [j[at + p x (t)]} exp [j6 x (l)] (115) 

and 

e,(r,<) = tf| exp [j[ a t + n (t)]\ exp Li«,(r)]. (110) 

B.I Amplitude Modulators 

When the waves given by (115) and (116) are incident on an analyzer 
whose axes are inclined at an angle 6 to e ± , the transmitted wave is 

e t (t,t) - c ± (r,t) cos 6 + e n (r,t) sin 6. (117) 

The total transmitted power is proportional to 



J | e t {r,t) \ 2 da, 



the integration extending over the beam cross section, and contains the 
ac terms or information carrying output terms 

P(t) | no <x E L Ei cos 6 sin cos <p(t) I cos 5(r) da 

! ■ 1 (U8) 

— sin <p(t) I sin S(r) da 
where ^(0 = tp x (t) - n (t) and 8(r) = 5 x (r) - 5„(r). 

B.2 Heterodyne Detection 

In this system we beat in a square-law mixer a phase-modulated sig- 
nal with a temporally and spatially coherent local oscillator, obtaining 
an output proportional to the product of the local oscillator and signal 
amplitudes. For a local oscillator wave of the form 

Cl . o (r,0 = Eio exp [ji^ot + *)] (119) 

the wave amplitude of the output difference frequency is given by 

e d (T,t) cc E ± , |,tf LO cos [(w- u>x.o)t- + + <p x , l (t) +« Xl ||(r)]. (120) 
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This represents the wave amplitude produced at each point on the sur- 
face of the mixer. What is of significance beyond this point is the aver- 
age amplitude (e d (t)), since this determines the response of the IF sys- 
tem.* Averaging over the mixer cross section A m gives 



(edit)) = E ± , i|J5 LO <cos[(a) - w LO )f - ^ + <p x , | (01 

X -r- / cos5 Xl I, (r) da - sin [(w - « LO )i — $ (121) 

A,„ Js 

1 f 1 

+ <P±. i(01 j- J sin p x , ,(r) do . 

B.3 Phase Modulator and Optical Discriminator 

To be specific, we consider an elementary discriminator formed by 
operating on the side of the resonance of a Fabry-Perot etalon. If we 
operate at the half-transmission point, the output wave is to first order 
given by 

«o(r,0 = |[1 + M( WlnBt - u)]e ± . ,(r,0 (122) 

where M is the slope of the transmission curve at co, and a>i n8t is the 
instantaneous frequency 



(123) 



toinat = -jr. [at + </>x , ii (0 + 5 x . n (r)] 

, °V± , || (Q 

= CO -+■ • 

^ dt 

The total transmitted power is proportional to 

j |eo(r,0 \ 2 da, 

and for small frequency deviations the output power is therefore given 
by 

P(t) -[l+M^i^J^.o^lexpb^.^r)]! 2 ^. (124) 



♦Regardless of the form of detector (photodiode, 39 - 40 photoconductor 41 or 
photoemitter 42 - 43 ) the optical mixing occurs over a finite area large compared with 
an optical wavelength, but small compared to the different frequency wave- 
length. Thus we consider that the total difference frequency output of the mixer 
is a sum of the independent contributions from each elemental portion of the ac- 
tive area, i.e., it is proportional to the value of ea averaged over the photo surface. 
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It should be borne in mind that (124) is strictly valid only when the 
discriminator is illuminated with perfectly collimated light. Strain has 
the additional complicating effect of producing a spread in the beam in- 
cident on the discriminator. Such a distribution of incident angles is 
equivalent to a broadening of the discriminator frequency response, and 
hence reduces the FM-to-AM conversion efficiency. 

appendix c 

The coupling structure analyzed in this appendix is the Fabry-Perot 
resonator shown schematically in Fig. 10(a). The dimensions of the 
resonator L c (length) and l c (spacing) are both assumed to be much 
larger than the optical wavelength. The structure can therefore be 
analyzed in terms of plane waves bouncing back and forth between 
plates. Since the length of the coupler is finite, we assume that energy 
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Fig. 10 — Fabry-Perot coupler for optically matching to the modulator when 
the aperture exceeds the reflector spacing: (a) geometrical arrangement; (b) equiv- 
alent circuit. 
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leaks off to the right in the y direction with group velocity v = (c/n ) 
sin 0, . Neither diffraction losses nor losses within the electro-optic me- 
dium are included in this treatment. 

We take as a model for the Fabry- Perot resonator a resonant trans- 
mission line system with equivalent circuit shown in Fig. 10(b). The 
resonator is represented by a transmission line which is an integral num- 
ber of half wavelengths long in the transverse or x direction. The "walk- 
off" losses are represented by assigning a small resistance, Z L = rZ 2 , to 
the terminating short circuit. Losses can be introduced in a variety of 
ways, but for high-Q resonators the results are independent of the way 
the losses are introduced; for the present treatment the representation 
chosen is particularly convenient. The finite reflectivity dielectric multi- 
layers are represented by an ideal transformer. Strictly speaking, the 
transformer turns ratio n will depend upon the angle of incidence and 
details of the multilayer construction. Experimentally, n is deduced by 
measuring the reflectivity of the multilayers when the resonator is 
terminated so as to eliminate reflections from the back surface (i.e., 
the perfect reflector). Finally, the characteristic impedances Z x and Z 2 
refer, respectively, to the wave impedances of free space and the electro- 
optic material in the x direction. They depend upon the index of re- 
fraction, angle of incidence, and direction of polarization of the optical 
signal. In order to illustrate the design principles for such a coupling 
system, we note that for waves polarized perpendicular to the plane of 
incidence in Fig. 10(a), Z x and Z 2 are given by 9 

Z, = t? sec 0/ (125) 

and 

Zt« -3-880 ft (126) 

Wo 

where tj = 377 ohms is the intrinsic impedance of free space and n is 
the refractive index of the modulator material. 

The objective of the analysis which is to follow is basically to calcu- 
late the dimensions of the Fabry-Perot coupler and the required re- 
flectivity for 100 per cent transmission of light, in terms of known pa- 
rameters. We begin by observing that for a given aperture D b the coupler 
length L c is given either by 

L c = D b (l - n 2 sin 2 0,-)"* (127) 

where use has been made of Snell's law in relating 0, to 0/, or by 

L c = ZVcos0,-. (128) 
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Two equations are required to specify L c , since it is possible for the 
angle of incidence required to achieve velocity synchronism to exceed 
the maximum achievable incident angle 9, determined by the limit 
angle relation 

sin 0, = l/n . 

It is found in fact that in order to achieve velocity matching the ex- 
ternal angle of incidence must satisfy the relation 

/ Wo 2 

and for the case n 2 > VV, (e.g., CuCl) we cannot employ the arrange- 
ment of Fig. 10(a) for which (127) applies. In that case it is necessary 
to use a more complex coupler such as the one given in Fig. 11 and use 
(128) in place of (127) to calculate L c . The spacing /„ can next be found 
from the resonance condition fi x l e = mir, which yields 



lc = 



wiAo 



for m = 0, 1, 2, 



(129) 



2/io cos 6i 

The calculation of the requisite reflectivity for 100 per cent light trans- 
mission into the modulator can be found in the following manner. We 
can couple all of the incident energy into the modulator by choosing 
the turns ratio n (i.e., the reflectivity of the multilayers) to obtain 
critical coupling, in which case the loaded Q, Q L , becomes half of the 
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unloaded Q, Q , of the resonator. Since the resonator is an integral 
number of half wavelengths long, the condition for critical coupling is 
equivalent to choosing n to satisfy the impedance matching relation 

n 2 = r§- 2 . (130) 

Disregarding for the moment the determination of r, which, as will be 
shown subsequently, can be expressed in terms of the unloaded resona- 
tor Q, we find that the multilayer reflectivity, or square of the input 
reflection coefficient when transmission line Z 2 is terminated [see Fig. 
10(b)], must be chosen so that 

Re = (Lzl\. (im,) 



\1 + r) ' 



Equation (130) has been used in obtaining this result. 

There only remains now the problem of determining the dimensionless 
constant r. The unloaded Q of the resonator can be found by observing 
that energy must leak off in the ^/-direction [Fig. 10(a)] with group 
velocity v g = (c/n ) sin 0< . The power flow to the right then is 

where W,/L c is the optical energy stored per unit length in the Fabry- 
Perot, and so the resonator Q(Q = o>W t /P) is simply 

Q = 2*- ^- e n„ csc 0,-. (132) 

Ao 

On the other hand, for the equivalent circuit of Fig. 10(b) the unloaded 
Q is given by 

Qo = mtr/2r. (133) 

The constant r can now be found by equating (132) and (133). Then 
by using (127) or (128) and (129) to eliminate L c and the integer m, 
we finally arrive at the formula 

(Vl - n 2 sin 2 0,- for n 2 < V^ 
r - 1 £ S in B< cos 6 A (134) 

2 Dh (cos Oi for n 2 > V? . 

It should be noted that the final results drawn from this analysis are 
valid only when the light is polarized normal to the plane of incidence. 
The method of analysis outlined applies quite generally, however, and 
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can be used for any plane of polarization provided that the character- 
istic impedances Zi and Zz can be determined. For design purposes the 
expressions of particular significance are (127), (128) and (129), which 
fix the dimensions of the coupler, and (131) and (134), which give the 
necessary reflectivity for transferring 100 per cent of the light into the 
modulator (neglecting, of course, diffraction losses). 

One final comment is in order concerning the output coupler. It is 
evident almost by inspection that the output coupler could be of pre- 
cisely the same form as the input coupler but would transmit only half 
of the light. As an alternative, we suggest using a scheme whereby the 
beam is simply allowed to diverge at the output of the modulator. If the 
beam divergence as determined by (49) is not excessive, then in principle 
all of the light could be collected by a lens, thus making the output 
coupler 100 per cent efficient. 
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